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Abstract - The limit n oo of the classical 0(n) 4> 4 model on a 3d film with free surfaces 
is studied. Its exact solution involves a selfconsistent Id Schrodinger equation, which is solved 
numerically for a partially discretized as well as for a fully discrete lattice model. Extremely 
precise results are obtained for the scaled Casimir force at all temperatures. Obtained via a single 
framework, they exhibit all relevant qualitative features of the thermodynamic Casimir force 
known from wetting experiments on 4 He and Monte Carlo simulations, including a pronounced 
minimum below the bulk critical point. 
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A celebrated example of fluctuation-induced forces is 
the Casimir force between two metallic, grounded plates 
in vacuum Such forces caused by the confinement 

of quantum electrodynamics (QED) vacuum fluctuations 
of the electromagnetic fields are expected to have con- 
siderable technological relevance. This has made them 
the focus of much ongoing research activity. During the 
past two decades, it has become increasingly clear that a 
wealth of similarly interesting classical analogs of such ef- 
fective forces, induced by thermal rather than quantum 
fluctuations, exist [3]j^] Two important classes of such 
"thermodynamic Casimir forces']^] are forces induced by 
fluctuations in nearly (multi) critical media between im- 
mersed macroscopic bodies or boundaries, and forces due 
to confined Goldstone modes [6]. Clear experimental ev- 
idence for the existence of such thermodynamic Casimir 



1 For a review of the Casimir effect in QED and an extensive lists 
of references, see |2| 

2 For reviews of" the thermodynamic Casimir effect and extensive 
lists of references, see [lj 

3 Following established conventions we use the term "thermody- 
namic Casimir forces" for forces induced by thermal fluctuations, in 
particular, also for near-critical Casimir forces, reserving the name 
critical Casimir forces to those where the medium is at a critical 
point. This topic must not be confused with those of thermal effects 
on QED Casimir forces and thermal Casimir-Polder forces, which are 
less universal since material properties of the media and confining 
objects matter; see, e.g., [5] 



forces was provided first indirectly by measurements of the 
thinning of 4 He wetting films at the A-point as the temper- 
ature T is lowered below the bulk critical temperature T c 
7|. Subsequently, direct measurements of the thermody- 
namic Casimir force on colloidal particles in binary liquids 
near the consolute point could be achieved [8j. 

Despite obvious analogies, crucial qualitative differences 
between thermodynamic and QED Casimir forces exist. 
First, the latter usually can be studied in terms of effec- 
tive free field theories in confined geometries where the 
interaction of the electromagnetic field with the mate- 
rial boundaries is taken into account through boundary 
conditions. By contrast, investigations of thermodynamic 
Casimir forces at (multi) critical points necessarily involve 
interacting field theories. Second, whereas electromag- 
netic fields average to zero in the ground state, the thermal 
averages (<p(x)) of fluctuating densities 4>(x) (order pa- 
rameters) associated with thermodynamic Casimir forces 
do not necessarily vanish. Nonzero profiles (4>( x )) can oc- 
cur at all temperatures T when the symmetry (p — ► — 4> i g 
explicitly broken (as it generically is for fluids and binary 
fluid mixtures in contact with walls), and at low temper- 
atures if this symmetry is spontaneously broken. When 
a nonzero profile exists, it will respond to a change of 
the separation between boundaries and hence cause an 
effective force even in the absence of fluctuations. The 
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presence of such nonfluctuating background contributions 
to thermodynamic Casimir forces substantially spoils the 
analogy with the QED case. Obviously, thermodynamic 
Casimir forces in 4 He films provide much better analogs 
than those in binary mixed fluids because neither a spon- 
taneous breakdown of the U(l) symmetry is possible for 
finite thickness L of the film, nor an explicit breakdown of 
the symmetry through the boundary planes. 

It is therefore very unfortunate that no viable the- 
ory other than sophisticated Monte Carlo simulations 
of 



lattice XY models [9-11 has emerged which is ca- 
pable to yield, within a single framework, all experi- 
mentally observed [7] relevant qualitative features of the 
reduced Casimir force pFciT^L), where /3 = 1/ksT. 
For a d-dimensional film, the latter takes the scaling 
form f3Fc{T,L) ~ L~ d {)(x) with scaling variable x = 
t(L/^+) 1/u on large length scales^ where t = T/T c - 1 
while v and £ + are the critical exponent and t > am- 
plitude of the bulk correlation length, respectively. The 
observed features of the scaling function i9(x) are: (i) 
d(x) < for all x, (ii) a relatively small critical value 
0(0), (iii) a smooth minimum at x m i„ < 0, and (iv) 
a nonvanishing T — > limit oo). A theory based 
on renormalization-group (RG) improved Landau theory 
12 13 , though giving some insight, suffers from severe 
deficiencies: it erroneously predicts (dl) an ordered low- 
temperature phase for d — 3 and L < oo, (d2) a much to 
deep minimum of $ whose derivative has a (d3) jump dis- 
continuity there, and (d4) a vanishing T — > limit oo). 

In this Letter, we report exact results for the Casimir 
force of the 0(n) 4> 4 model on a film ]R 2 x [0, L] with free 
surfaces at z = and z = L in the limit n — ¥ oo. Its 
n = 2 analog describes 4 He fluid films near the A-point. 
The model has, for general n > 2, an ordered low-T bulk 
(L = oo) phase. For finite L, long-range order is restricted 
to T = since low-energy (spin wave) excitations destroy 
long-range order at any temperature T > 0. Further, the 
presence of confined Goldstone modes at T = implies a 
nonzero T — > limit of the Casimir force. Owing to the 
breakdown of translation invariance along the z-direction 
perpendicular to the surfaces, the exact n — > oo solution 
does not correspond to a mean spherical model with global 
constraint [M], but involves a Id Schrodinger equation 
with a selfconsistent potential V(z). The exact scaling 
function $(x) for n — oo can be expressed in terms of the 
eigenvalues and eigenfunctions of this equation. Deter- 
mining them by numerical means, we managed to get the 
extremely precise results displayed in Fig. [T] Determined 
within a single theoretical framework, these exhibit all fea- 
tures (i)-(iv) mentioned above but none of the deficiencies 
(dl)-(d4). 

Our best estimate for the Casimir amplitude is 

A c = tf(O) /2 = -0.01077340685024782(1). (1) 



4 Throughout this work, the symbol ~ 
cally equal" in the respective limit, e.g., 
li miwoo f(L)/g(L) = 1. 
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Fig. 1: (Color online) Scaling function i?(a;) of the critical 
Casimir force for L = {65,97,129,193,257} (model A, solid 
lines) and L = {97, 129, 193, 257} (model B, dash-dotted lines, 
for model definition see text). The dashed horizontal line in- 
dicates the Goldstone value #(— oo) = — f(3)/87T, and the dot- 
ted curve represents the limiting behavior &(x) ~ $(— oo) + 
(cln|x| + d)/x, with c = 0.07856(5) and d = 0.3673(5). The 
effective thickness L e ff accounts for leading scaling corrections 
(see text). Note that deviations from the excellent data col- 
lapse are only visible under magnification. 



The function 0(x) has a pronounced minimum i9 m i n = 
-0.1268565841360(1) at a; min = -4.55702477008(1). In 
the low-temperature limit x — > — oo, it approaches the 
value -C(3)/8tt = -0.0478283245 ... [15] of a massless 
free theory with Neumann boundary conditions^] Using 
an appropriate low-T model — a nonlinear a model on a 
film with Neumann boundary conditions [16) — we can 
show that the asymptotic behavior is of the form d(x) 
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tf(-oo) ~ 

Both the form of these logarithmic anomalies and the 
existence of a minimum at a: m i n < are intimately related 
to the breaking of translational invariance across the slab. 
To appreciate this, one should note the following. As may 
be gleaned from [18] and will be shown below, the exact 
n — > oo solution reduces to a constrained Gaussian model 
with an effective 4> 2 /2 term whose interaction constant 
f— V(z) involves a selfconsistent potential V(z). Accord- 
ing to [18], the t = analog of V(z) for the semi-infinite 

-l/4z 2 
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This algebraic be- 



case L = oo is given by 
havior carries over to the near-boundary behavior of V{z) 
for t < on scales small compared to both the Joseph- 
son coherence length oc \t\~ 1 ' and L 17 . By contrast, the 



and L 

counterpart of V(z) for periodic boundary conditions is 
independent of z, an enormous simplification which en- 



5 Arguments similar to those used in [16| show that the low-T 
limit is described by a nonlinear a model on a film with Neumann 
boundary conditions. 
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ables one to determine the correspon ding scaling func- 
tion $ pbc (a;) in closed analytical form [19]. It decreases 
quickly and monotonically from zero to its Goldstone value 
t? pbc (— oo) = — e(3)/V, has no minimum at finite x < 0, 
and approaches the x — > — oo limit ~ \x\ e - '^'. 

We next turn to an outline of the essentials of our cal- 
culations. The model is described by the Hamiltonian 
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with d = 3, where tfr = (4> a ) is an rt-component field and 



v e 



pd-l 



denotes the lateral coordinates. The boundary 



terms entail the boundary conditions (d z — ci)0| z= o and 

(9* +C 2 )0U=i. 

The n — > oo limit can be derived by standard means. 
Upon making a Hubbard-Stratonovich transformation, the 
partition function Z = J D[<p] e - ^ can be written as 



Z oc /£>[</>] e" 



!/d d - I s,/ i d*[<«f+iV>-V 2 )<*.- 



] , (3) 



where the Laplacian V 2 is subject to the above-mentioned 
boundary conditions. In the limit n — > oo, the ^-integral 
can be evaluated by saddle-point integration. Writing the 
saddle point as iipo(z) = V(z) — f , the reduced free energy 
per area A = J d d ^ 1 y^ Jl = —(An)^ 1 \nZ, becomes 



h = \l I 



^^{z\ln(j?-dl + V)\ Z ) 



T-V(z) 



AO) 
Jl ' 



(4) 



where is a trivial background term which does not 
matter henceforth. 

The stationarity condition 5Jl/5V(z) = yields 

f v(z) _ 9 f d^p ^ \M*)\ 2 (5) 



where e„ and ¥>i/(z) = are the eigensolutions of 
[-& Z +V(z)]<p v (z) = e v <p v {z). 



(6) 



The p-integrals on the right-hand sides of Eqs. Q and 
([5| are ultraviolet (UV) divergent at d = 3. To make the 
model ([2]) well defined and suitable for numerical calcula- 
tions, we must regularize these divergences. We study two 
distinct regularized versions of model In the first (A), 
only the z-coordinate is discretized, and the p-integrals 
are regularized dimensionally. In the second (B), a fully 
discrete lattice model is investigated. 

Model A — The system consists of L layers located at 
z = 1, • • • ,L, where we replace the operator d\ in Eq. ([6| 
by its discrete analog, the LxL matrix D 2 = (— 26 ZiZ ' + 



5\ z - z i\ } i). Rather than including analogs of the boundary 
terms oc 6j, we impose the Dirichlet boundary conditions 
</>| z=0 = 4>\ z= l+i — 0. The thickness change L — » L + 1 
is accounted for in the numerical analysis by introduc- 
ing an effective thickness L e g (see below). The Hamilto- 
nian of Eq. (l6| becomes the matrix H = — D 2 + V with 
V = diag(Vj)]_] The dimensionally regularized p-integral 
in Eq. ^ is straightforward. It produces a simple pole at 
d = 3, which gets absorbed in the bulk critical value f c . To 
see this, we subtract from Eq. ^ its bulk analog, making 
the appropriate replacements J2v ■ ■ ■ Wvz^ 2, JZ an d 
t v — > e(k) = 7'b + 4sin 2 (fc/2), where r b is the inverse bulk 
susceptibility. Writing f — t c + t and noting that r b = 
at t c , we can set d = 3 to obtain 



24ng-\T-V z ) = (z\]nH\z) 



(7a) 



To eliminate the UV singularities of the bulk free energy 
density / b — lirnL_ > . 00 /l/L we subtract from it its Taylor 
expansion to first order in r. This gives an UV finite 
renormalized fb(r,g) which at d — 3 becomes 



/b(r, g) = ^ V^b(4 + r b ) - 2+ ^ h arsinh yJr h /A 



-f (r-r b ) 2 , 



(7b) 



where r b is the solution to r b = t — j|- arsinh y/rb/4 
or depending on whether r > or < 0. The chosen 
bulk counterterms also absorb the UV singularities of /l. 
Their contributions cancel in the excess free energy / cx = 
h 



Lf b , yielding the UV finite d — 3 result 

fex(T,g,L) 

- Lf h (T,g) 
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i-tr[H(l-lnH)]- A tr[(r - V) 2 ] 



(7c) 



Solving Eq. (7a) numerically first for t = 0, we com- 
pute /cx(0,.9, L) to determine the surface free energy 
/ s (0, g) = /ex(0, g, oo)/2 and the Casimir amplitude Ac = 
limi_ i . 00 L 2 [f ex (0,g, L) — 2/ s (0, g)}. To extract precise val- 
ues from the data, knowledge about corrections to scaling 
is important. Clearly, the usual Wegner corrections gov- 
erned by the n — ¥ 00 exponent uj = A—d must be expected. 
Further, deviations from asymptotic Dirichlet boundary 
conditions are known to be described by irrelevant sur- 
face scaling fields Xj oc 1/ 6j ( "extrapolation lengths" ) that 



scale naively 21 p22 , so that their correction-to-scaling ex- 
ponent is u)\ = 1. Noting the degeneracy uj = ui\ = 1 
at d = 3, standard RG considerations can be used to 
show that leading corrections to scaling oc L^ 1 InL along 
with those oc L^ 1 should occur for the effective amplitude 
A c (g,L) = L 2 [f CK (0,g,L) - 2f s (0,g)]. As can be seen 



from the results for Ac depicted in Fig. [ZJ convergence is 
very poor when g is small, e.g., for g = 1. Relying on the 
results for Ac(l,£ < 64), one could easily infer an incor- 
rect value of —0.026(1) for Ac, more than twice as large 



S A preliminary analysis of the critical case was made in [20| 
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Fig. 2: (Color online) Effective Casimir amplitude Ac(g, L) for 
different values of g = {1, 2,4,..., 512, oo}. The solid lines are 
fits including corrections oc L _1 lnL, which vanish as l/g for 
g — > oo (see text). The dotted lines are guides to the eyes. 



as the correct asymptotic one Ac(l, oo) = —0.0108(1), see 
Fig. [2| Convergence is far better for g 3> 1 , as the min- 
imum in Ac is at L mm w 80/ g. Fast convergence and 
high precision can be achieved by solving Eqs. Q with 
r = and g set to its fixed-point valueQg* = oo, where 
the corrections oc L^ 1 InL vanish. 

To extend this g — > oo analysis to t ^ 0, we note 
that the critical exponent v = 1 and absorb the am- 
plitude £+(<?) = <?/247r of the bulk correlation length 



A+) 



-1/2 



£+r v for T > T c in the temperature 



variable, defining t = r/£+(g) so that the scaling variable 
becomes x = tL. In this limit, Eqs. ^ reduce to 



t= (zIlnHl 



z 



fh(t) 



1 J sinht -t for i > 0, 
4tt 1 for f < 0, 



(8a) 



(8b) 



/«(*,£) 



l 



i/b(t). (8c) 



-tr[H(l + t-lnH)]-g 

From Eqs. ([8]) we calculate the Casimir force fiFcit, L) = 
-dU(t,L)/dL « -[/„(t,L+ 1) - /«(*,£- l)]/2. To 
achieve the excellent data collapse shown in Fig. [l] it turns 
out to be sufficient to write -d(x) ~ Lg ff (3Fc(t, L), intro- 
ducing an effective thickness L c ft = L + SL, as proposed 
in 24 and substantiated by field theory [21]. 

Inspection of Eq. (8a) reveals that the 
eigenvalue cqL 2 is positive for all x and vanishes 



scaled lowest 



xe 



7 Recall that the bulk analog of the integral in Eq. (jSJ contains a 
rcgularization-depcndcnt term, which reads — a(d)r^ in the notation 
of Ref. [23] . Its coefficient a(d) may be identified as l/g* if a(d) > 0. 
For our regularization method A, a(d) = as n — > oo, so that 
g* = co. For method B, a(d) < 0. Hence no infrared-stable g* > 
exists and corrections to scaling cannot be suppressed in this manner; 
cf. Ref. [23] Sect. 2.4]. 



as x — > —oo. That is, for finite L < oo, the system remains 
paramagnetic whenever T > 0, due to the nonperturba- 
tive generation of a mass. Furthermore, the remaining 
eigenvalues e v> o approach the Neumann values. 

The numerical calculations are performed with 33 dig- 
its precision. This yields about 30 significant digits 
in / cx . For the effective thickness the form L e g = 
L + SL + YliLi biL~ l is chosen. The estimates of Ac 
and SL are then determined by analyzing f ex (0,L) for 
L = 1600, 1800, . . . , 3800 and L = 4096 with the ansatz 
f„(0,L) = 2/ s (0) + A c Lj. Our final results, Eq. 
and SL = 0.7255032704723(3), are obtained using m = 5. 
As benchmark for the errors, the variations of the esti- 
mates resulting from analogous analyses with m = 4 and 
different choices of thicknesses L are used. 

Model B — A simple cubic lattice model of sites x = 
(y,z) G Z d with Xi = l,...,Ni is considered whose 
Hamiltonian follows from Eq. ([2| through the replace- 
ments / d^y J Q L dz -> £ x and (V0) 2 -> £<=i |0(x) - 
4>(x + e.j ) | 2 , where are unit vectors along the principal 
axes. Along the z-direction, Dirichlet boundary condi- 
tions are again imposed; along all y^-directions, periodic 
boundary conditions <f> y = <fi y+ N iei are chosen and the 
limits Ni <c i — > oo taken at fixed L = N^. Just as in model 
A, the operator d 2 in Eqs. Q and ^ is replaced by its 
discrete analog D 2 . Owing to the discreteness of the lat- 
tice, the momentum integrations are restricted to \pi\ < tt. 
Further, to account for the modified dispersion relation, 
4 sm2 (Pi 1 2) must be substituted for p 2 . Thus the p- 

integral of each series coefficient of the spectral sum in 
Eq. ([5]) becomes an analytically computable e„-dependent 
two-dimensional Watson integral. Its three-dimensional 
bulk analog can also be determined analytically 25 , along 
with their antiderivatives one encounters in the analogous 
momentum integrals of the right-hand side of Eq. Q . 

For model B no fixed-point value g* > exists to 
which g could be set to eliminate leading scaling cor- 
rections 23. Thus, scaling corrections oc L _1 lnL re- 
main for any value of g. However, the amplitudes of 
these corrections become minimal in the limit g — ¥ 00. 
We use L c ff = L + a lnL + SL + J2iLi( a i^ n ^ + h)L~ i 
with m — 3 to analyse / cx (0,i), which we computed for 
various thicknesses up to L = 4096. We reproduce the 
value ([!]) for the Casimir amplitude to 13 significant dig- 
its — a striking confirmation of universality. Furthermore, 
we get a = -0.123903101(1) and SL = 0.81422072(1). 
Analyzing our data also at the minimum of the thermo- 
dynamic Casimir force, we obtain results for x m i n and 
$min that are fully consistent with those obtained from 
model A. The value of ao is consistent with the one ob- 
tained above from the analysis of f cx (0,L). However, we 
get SL = 1.1979(1) and SL = 0.8924(1) from the loca- 
tion and the value of the minimum of the thermodynamic 
Casimir force, respectively. The plot of $ in Fig. [T] uses 
L e g = L — 0.1239031 lnL + 1, which results in a good data 
collapse. However, we should keep in mind that this way 
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corrections oc L _1 cannot be eliminated completely for the 
whole range of the scaling argument x. 

Finally, let us comment on the stability of our results 
with respect to a change of boundary conditions, as in- 
vestigated recently in 28 . For both models A and B, the 
nearest-neighbor (NN) bonds @Jj in the two boundary 
layers j = 1 and j — L were chosen to agree with the NN 
bonds fij in all other layers. Allowing for boundary cou- 
pling^] Ji and Jl ^ J, one can change the corresponding 
diagonal elements of D 2 . However, this change is com- 
pensated by a corresponding change of the potential V. 
As expected for d = 3, the self consistent solutions remain 
asymptotically the same. 

In summary, considering the 0{n) vector model on a 
slab with free boundary conditions, we expressed its uni- 
versal scaling function for the Casimir force in the limit 
n — > oo exactly in terms of the eigensystem of the re- 
sulting selfconsistent Id Schrodinger equation, which we 
then solved by numerical means using two qualitatively 
distinct regularization methods. We obtained consistent 
results that agree to many digits. They exhibit all qualita- 
tive features (i)-(iv) one expects to hold for general n > 2, 
in particular, nontrivial crossovers from three-dimensional 
critical to two-dimensional pseudo-critical behavior and to 
the low-tcmpcrature behavior dictated by confined Gold- 
stone modesj^] Besides being interesting in their own right, 
these results could pave the way to successful approxi- 
mate analytical treatments of such challenging problems. 
We expect them to play a role as fruitful as large-n solu- 
tions have done in the theory of quantum critical phenom- 
ena [29] . Furthermore, our analysis can be extended in a 
straightforward fashion to include bulk and surface mag- 
netic fields and to study appropriate quantum versions of 
the model. 
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